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Abstract. We present a reconstruction algorithm for recovering both "magnetic- 



hard" and "magnetic-soft" obstacles in a background domain with known 
isotropic medium from the boundary impedance map. We use in our algo- 
rithm complex geometric optics solutions constructed for Maxwell's equation. 



1. Introduction 



In this paper we study an inverse boundary value problem for Maxwell's equa- 
tion. Let C I 3 be a bounded domain with smooth boundary, filled with 
isotropic electromagnetic medium, characterized by three parameters: the per- 
mittivity e(x), conductivity cr(x) and permeability jx{x). A "magnetic-hard" ob- 
stacle is a subset D of £1, with smooth boundary, such that the electric-magnetic 
field (E, H) satisfies the following BVP for Maxwell's equation 



where v is the unit outer normal vector to the boundary dflUdD. The boundary 
condition v AH\qd = physically appears when an active object is presented. 
Another type of obstacle, what we call a "magnetic-soft" obstacle, is a subdomain 
D such that the tangential component of the electric field v A E vanishes on the 
interface dD, when a passive object is presented. Define the impedance map by 
taking the tangential component of the electric field v A E|gn to the tangential 
component of the magnetic field v A Hgc^. Then our purpose is to retrieve infor- 
mation of the shape of D from the impedance map. 

The well-known Calderon's problem [2] is to determine the conductivity of a 
medium by making voltage and current measurements of the boundary. The 
information is encoded in the Dirichlet-to-Neumann map for the conductivity 
equation V • (7V?/) = 0. In [16], Sylvester and Uhlmann constructed complex 
geometric optics (CGO) solutions for Schrddinger operator A — q and proved 
the uniqueness of C 2 isotropic conductivity in dimensions n > 3. Further de- 
velopments including improved regularity assumption, 2D problems and partial 
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data problems, were obtained based on this idea. See [18] for a survey of recent 
developments. 

Another application of CGO solutions, the enclosure method was first intro- 
duced by Ikehata [5, 6] to identify obstacles, cavities and inclusions embedded 
in conductive or acoustic medium. Geometrically, using the property of CGO 
solutions that decay on one side and grow on the other side of a hyperplane, 
one can enclosing obstacles by those hyperplanes. This idea was generalized to 
identify non-convex obstacles by Ide et al. [7] for isotropic conductivity equa- 
tions (conductive medium) and by Nakamura and Yoshida [11] for Helmholtz 
equations (acoustic medium), by utilizing the so-called complex spherical waves 
(CSW), namely, CGO solutions with nonlinear Carleman limiting weights. In 
[19], Uhlmann and Wang constructed generalized complex geometric optics so- 
lutions for several systems with Laplacian as the leading order term, e.g., the 
isotropic elasticity system, and implemented them to reconstruct inclusions. 

As for Maxwell's equation, [15], [12] and [13] answered the uniqueness question 
for parameters with suitable regularity from the impedance map in a domain Q. 
In [13], the Maxwell's operator was reduced into a matrix Schrodinger operator 
and vector CGO solutions were constructed to recover electromagnetic parame- 
ters. 

To address the inverse problem of determining an electromagnetic obstacle, we 
observe that solutions of a non-dissipative Maxwell's equation (<r = 0) share sim- 
ilar asymptotical behavior (a key equality in Lemma 4.2) to those of Helmholtz 
equations (a key inequality in Lemma 4.1 in [5]). Therefore, with CGO solutions 
at hand, the enclosure method is applicable: one can define an indicator function 
Ip(r,t) for each direction p £ S 2 ; by adjusting t, the hyperplane moves along p; 
for each p and t, the asymptotical behavior of I p (t, t) as r 3> 1 produces the sup- 
port function of the convex hull of D. However, unlike the Schrodinger equation, 
the CGO solution for Maxwell's equation doesn't behave as small perturbation 
(w.r.t. r) of Calderon's solutions. We overcome this by carefully choosing rela- 
tively large "incoming" fields (w.r.t. r) compared to the perturbation. 

The rest of the paper is organized as following. In Section 2, we formulate 
the forward problem for a "magnetic-hard" obstacle and define rigorously the 
impedance map. Then we construct the CGO solution of interests in Section 3. 
The main reconstruction algorithm for "magnetic-hard" obstacles is introduced 
and proved in Section 4. Finally, we remark in Section 5 that the scheme also 
applies to "magnetic-soft" obstacles. Through the whole paper, well-posedness 
of a mixed boundary value problem for Maxwell's equation plays an important 
role. Hence, for completeness, we include a proof in Appendix A. 
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2. Direct Problems and the main result 

Let f2 be a bounded domain in IR 3 with smooth boundary and its complement 
IR 3 \ f2 is connected. We consider electric permittivity e(x), conductivity cr(x) 
and magnetic permeability (J,(x) of the background medium as globally defined 
functions with following properties: there are positive constants £ m ,£M, fJ-m, fJ-M, 
a Mi £o and no such that for all x G 

£m < S(x) < £ M , /JL m < fx(x) < S M , < a(x) < G M 

and e — eo, a, ji — /io G Cq{£1). 

A "magnetic-hard" obstacle D (corresponding to the sound-hard obstacle for 
Helmholtz equations) is a subset of 0, such that $7 \ D is connected. Moreover, 
the electric field E and the magnetic field H in £l\D satisfy the following boundary 
value problem of the time-harmonic Maxwell's equation 

V A E = iufM, V A H = -iwyE in £1 \ D, 
uAE = fe TH^ v {dn) on dtt, ' L ' ' 

where 7 = e + i~, and the "magnetic-hard" boundary condition on the interface 
dD 

(MH)| Sfl = 0, (2.2) 

where v is the unit outer normal vector on dVL U dD. Through out this note, we 
assume the non-dissipative case a = 0. Then 7 = e is a real function. 

Notations. If F is a function space on <9S7, the subspace of all those / G F 3 
which are tangent to <90 is denoted by TF. For example, for u G (fP($S7)) 3 , 
we have decomposition u = ut + u u u, where the tangential component ut = 
— v /\{y Au) G TH s (d£l) and the normal component u u = u ■ v G H s (d£l). For a 
bounded domain Q in IR 3 , we denote 

THfy. w (dn) = {/ G TF s (dft) I Div(/) G H s (dQ)}, 

H Div^ = {u G (F fc (^)) 3 I Div(^ Au\ dn ) G F*- 1 / 2 ^)}, 
with norms 

ll/llr^cen) = \\f\\ 2 H«m) + H Div (/)ll^), 
Nffft. («l) = H u llff*(n) + H Div (^ /\u\ dn )\\ 2 Hk _ 1/2{m) , 

where Div is the surface divergence. There are natural inner products making 
them Hubert spaces (see [15]). In addition, we define the weighted L 2 space in 
IR 3 : 

L 2 S = {/ G Lf oc (R 3 ) : H/ll 2 -^ J(l + \x\ 2 ) s \f(x)\ 2 dx<^y 
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1/2 

Admissibility. It can be shown (see Appendix A.) that for / G THj^^(dQ) and 

1/2 

g G TH-Q-_^{dD), the boundary value problem for Maxwell's equation 

V A E = iujfjH, V A H = -iur/E in U\D, 
uAE\ an = f (2.3) 
v AH| 9D =5, 

has a unique solution (E, H) G iJ^. (0 \ D) x H^ v (£l \ D), except for a discrete 
set of magnetic resonance frequencies {uJ n }- It satisfies 



ll E H^ iv (^) + ll H H^ iv (^) ^ C ^/llr^ v («i) + ^th^d) ) ■ ( 2 - 4 ) 

Denote by (Eq, Ho) the solution of Maxwell's equation in the domain 0, without 
an obstacle, namely, the solution of 

VAEo = iw/^Ho, V A Ho = — iweEo, in O, 



Main result. Now we are in the position to define the impedance map for 
non-resonant frequencies, 

A D (uAE\ an ) = v AH| 9n , 

and it can be shown that 

a d : TH^an) - th^xi) 

is bounded. If u is a resonance frequency, one can replace the impedance map 
by the Cauchy data set 

C u = {{v A Elan, v A E|an) I (E, H) satisfies (2.1) and (2.2)} 

C TH^dn) x TH^(dn). 

Denote by A0 the impedance map for the domain without an obstacle. Then the 
main result of the presenting work is to show 

Theorem 2.1. For non-dissipative background medium (a = 0), there exists a 
reconstruction scheme for the obstacle D from the impedance map Ap. 

3. Complex geometric optics solutions 

In [13] , the Maxwell's operator was reduced to an 8 x 8 second order Schrodinger 
matrix operator by introducing the generalized Sommerfeld potential. A vec- 
tor CGO-solution was constructed for the Schrodinger operator, simplifying the 
proof in [12]. Similar techniques also appeared in [3] when dealing with the in- 
verse boundary value problems for Maxwell's equations with partial data. For 
completeness, we include the construction of the solution in this work (see [13] 



for more details) and provides a special choice of the "incoming" fields. 



Define the scalar fields <5 and ^ as 

i 



V-( 7 E), * = -V-(/xH). 



(3.1) 



Under appropriate assumptions on $ and ^F, Maxwell's equation is equivalent to 



V AE 



1, 



-V 



1 



1, 



1 



iuJn~H = 0, VAH + -V -$ +iw7E = 0. (3.2) 
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^ J A* 
Moreover, in this case, $ and \& vanish, leading to a solution of Maxwell's equa- 
tion. Let X = (ip,e,h,i;) T £ (V'f with 

e = 7 1/2 E, h = Ai 1/2 H, 



1 



1 



(3.3) 



P(iV) 



F=(fc-«)1 8 + 



\ 



D 



D 



-i 



/ 



Then (3.1) and (3.2) read 

(P(»V) - k + V)X = 0, in n 

where 

f V- \ 

V VA 
-VA V 

V o o v- o / 

/ / V- \_ 

V -VA 
VA V 

V o o v- o / 

are matrix operators and 

D = diag( Ai 1 / 2 ,7 1/2 l 3 ,A i 1/2 l3,7 1/2 ), « = «(7/i) 1/a , fc = "(W- 

An important property of this operator is that it allows to reduce Maxwell's 
equation to the Schrodinger matrix equation by noticing 

(P(iV) - fc + V)(P(iV) + fe - V T ) = -(A + k 2 )l 8 + Q, (3.4) 

where 

Q = VP(iV) - P{iV)V T + k(V + F T ) - VV T 

is a zeroth-order matrix multiplier. Hence, by writing an ansatz for X, one can 
define the generalized Sommerfeld potential Y 

X = (P(iV) + k — V T )Y. 

So it satisfies the Schrodinger equation 



(-A-k 2 + Q)Y 
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0. 



(3.5) 



The following CGO-solution is constructed using Faddeev's Kernel. Let ( £ C 3 
be a vector with ( • ( = k 2 . Suppose yo,C € C 8 is a constant vector with respect 
to x and bounded with respect to £. We refer e ?x '^yo,c as the "incoming" field. 
Then there exists a unique solution of (3.5) of the form 

Y c (x) = e ix <(y , c -v c (x)), 

where v^(x) G (L 2 +l ) 8 satisfying 

\h\y s+1 < c/ici 

for 5 G (—1,0). Moreover, one can show that S (-fP(O)) 8 for < s < 2, e.g., 
see [1], and 

\\v c (x)\\ H . {a) <C\{\-\ (3.6) 



Lemma 3.1 in [13] states that if we choose yo,c such that the first and the 
last components of (-P(C) — ^)z/o,c vanish, where -P(C) is the matrix obtained by 
replacing iV by C in f(£V), then for large |£|, -X^ provides the solution of the 
original Maxwell's equation. We proceed to provide a more specific choice of yo.c 
such that the CGO solution for Maxwell's equation has special properties. 



As in [13], choose 



where 



2/o,c = T^t(C • a,ka,kb,(- b) T , 



C = —irp + \f t 2 + k 2 p ± , 



with p, p 1 - G S 2 and p-p L = 0. r > is used to control the size of |C| = V2r 2 + k 2 . 
Then we obtain 



xox ■= (P(-0 + %o,< 



ICI 



/ \ 

-{(■a)C-K Ab + k 2 a 
kC A a - (C • b)C + k 2 b 
V 



satisfying the condition in Lemma 3.1 in [13]. 
Taking r — » oo, we have 



We choose o and b such that 



C * 6 = l, C-« = o. 
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This is satisfied, for example, by taking a G M 3 and b G C 3 satisfying alp, 
a 1 p 1 and b = (. Given these choices, It's easy to see that 

i] := (xo,c)2 - ► — A 6 = i/cp A p 1 - (~ 0(1)) as t — > cx), 

:= (3:0^)3 ~ 0{t) as t — > 00. 
Then Xq = (P(iV) + fc — y T )Y^ is written in the form 



X< = e T(x.p)+iVT2+k2x.p±( XQ ^ + 

where 

r<; = p (-()vc + P(W)v c - V T y ^ + kv ( - V T v c (3.7) 
satisfying for C > independent of £ 

Ikclli^(n) ^ c '- 

Summing up, we obtain the following 

Proposition 3.1. Let p, p 1 - G §> 2 mi/i p • p -1- = 0. Assume u is not a resonant 
frequency. Given 9, i] G C 3 as above, then for r > /ar^e enough, there exists a 
unique complex geometric optics solution (Eo,Ho) of Maxwell's equation 

V A E = iupHo V A H = —iLoeF, in 0, 

of the form 

E = £ -l/2 e T(x-p)+»v/r 2 +Px./( J) + R j 
H = ^/2 e T(x-p)+iVT 2 +k2x-p^ + 

Moreover, we have 

77 = 0(1), = 0(r) /orr» 1, 
and i? = (7^)2, Q = ( r c)3 are bounded in (L 2 (S7)) 3 /or r>l. 

For reconstruction, one needs to compute the boundary tangential CGO-fields 
(y A Eo|an, v A Ho|an). In [12, 13], by solving a boundary integral equation, one 
can recover them from the impedance map A0 if the background parameters are 
unknown. In our case with known medium, CGO-fields are known. 

4. Reconstruction Scheme 
Adding a parameter t > into the CGO-solution in Proposition 3.1, we use 

E = e -l/2 e r(x.p-t)+iVr 2 +fc 2 x.p^ (jj + R) 

to define an indicator function which physically measures the differences between 
the energies required to keep the same boundary CGO electric field for the domain 
fl with and without the obstacle D. 
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Definition 1. For p £ § 2 , r > and t > we define the indicator function 



where Eo is a CGO solution of Maxwell's equation given by (4.1). 

The enclosure method's aim is to recover the convex hull ch(D) of D by re- 
constructing the following support function. 

Definition 2. For p £ S 2 , we define the support function of D by 



Then, the reconstruction scheme in Theorem 2.1 is 

Theorem 4.1. We assume that the set {x € IR 3 | x • p = /i£>(/?)} n dD consists 
of one point and the Gaussian curvature of dD is not vanishing at that point. 
Then, we can recover ho(p) by 



Moreover, if D is strictly convex (the Gaussian curvature is everywhere positive), 
then we can reconstruct D. 

Remark 1. The proof of Theorem 4.1 mainly consists of showing the following 
limits: 



Remark 2. Gaussian curvature at a point on a surface is defined to be the product 
of two principal curvatures, which measure how the surface bends by different 
amounts in different directions at that point. A surface with positive Gaussian 
curvature at a point is locally convex. Note that the non-vanishing assumption 
on the Gaussian curvature in the theorem is not crucial since only finitely many 
directions p violate the condition. 

4.1. A key integral equality. To show the limits in Remark 1, we need the 
following equality for non-dissipative Maxwell's equation 

Lemma 4.2. Let a = 0. Assume (E,H) is a solution of 

V A E = iupH, V A H = -iooeE, in tt\D 

satisfying the boundary condition 




hn(p) '■= sup x ■ p. 



h D (p) = ini{t G R | lim I p (r, t) = 0}. 



lim ipfr, t) = 0, when t > hoip)', 
liminf/ p (r, h D (p)) = C > 0. 



isAH\ dD = and v A E| an = v A E |aa- 
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Then we have 



iu A E ) • (vAH-vA H ) A v 
Jdn 



dS 



/ /i~ X |V AE- V AE | 2 -cA|E-E | 2 dx 

Jf2\D 



+ / /i _1 |V AE | 2 -io 2 e\B \ 2 dx. 

ID 



(4.2) 



Proof: Denote 



I:=iu f (v A E ) • [{v A H - v A H ) A v]dS = iu f (u A E ) • (H - H )d5. 
Jan Jan 

First by integration by parts, we have 

/ A E) • (V A E - V A E ) - to 2 eE ■ (E - E )dx 

= -(/ -/ ) (1/ A /i _:L (V A E)) • (E - E )dS = 

by boundary conditions. Adding this to the following equality 

I = [ [y A E ) • (-iwH + iwHo)dS 
Jan 



-/z _1 (V A E ) • (V A E) + u/eE • Ecte 



+ / /i _1 |V AE | 2 - oo 2 e\E Q \ 2 dx + / A Eo) • (—iutVjdS 
Jn JdD 

with the last term vanishing due to the zero-boundary condition on the interface, 

/ (v A E ) • ( = kM)dS = [ (u A E ) • (—iuj(v A H) A v)dS = 0, 
JdD JdD 

we obtain (4.2). ■ 



4.2. Proof of Theorem 4.1. We proceed to show the first limit in Remark 1 

lim I p (T,t) = if t > h D (p) (4.3) 

T— +00 

by proposing an upper bound of the indicator function. 



Let E = E — Eo be the reflected solution in 0, \ D. It satisfies 

V A (/z -1 V AE) - s£ -1 V(V • eE) - u 2 eE = in Q \ D, 
v A E^ = 0, 

v A (//^V A E)\ 9D = -iuu A Ho|aD G TH^idD). 
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(4.4) 



By (A. 10) in Appendix A., we have 

\Mm(u\D) <C\W^ m 2 H -^8D) ^ c (ll v A E °\\h D ) + \\Vo\\h (D) ), (4.5) 
where the second inequality is valid since 

iu(u A H , F) dD = (/z^V A E , V A F) D - {iojeE , F) D for F G {H^D)) 3 . 
Therefore, by (4.2) and (4.5), we have 

I p (r, t) < C(||Eo||| 2(D) + || V A E f L2(D) ) < C(||Eo||| 2p) + \\H \\ 2 L2{D) ). (4.6) 
Plugging in the CGO-solution (4.1), we obtain the following estimates: 
||E ||| 2(D) < CeMMp)-*)^ + R\\ 2 3 ~ e 2 ^ D (p)-t) T » i, 
||H ||i 2(D) < Ce 2 ^Mp)-t)||0 + Q ||2 2 ^ )3 ^ T 2 e 2T(h D ( P )-t) T » i ( 4 - 7 ) 

Therefore, we obtain 

I p (T,t) < C T 2 e ^(^(p)-*) 
for t large enough, proving the first limit (4.3). 

To show the second limit 

liminf lJT,h D (p)) = C >0, (4.8) 
it suffices to prove the following two lemmas. 
Lemma 4.3. If t = ho{p) in CGO-solution (4.1), then 

liminf/" fj,- l N AE | 2 dx = C, 
T ^°° Jd 

with some constant C > 0. 

Lemma 4.4. If t = ho{p), then there exists a positive number c such that 
^j n ^e\E-E \ 2 dx + f D e\E \ 2 dx^ 



LO 2 



J D fi~ l \V /\E \ 2 dx 

for t large enough. 



< c < 1, (4.9) 



Proof of Lemma 4.3: This is obtained by noticing, in Proposition 3.1, that 
the first order growth of the constant vector 9 in Ho with respect to r. Then the 
left hand side integral 

J pT l \V A E | 2 dx > CHHoll^js >C J T 2 e 2T(x - p - hD{p)) dx >C for r > 1. 

To show the last inequality, we denote by xq the point in {x £ M 3 | x ■ p = 
ho(p)} H 3D. It's not hard to see that there exist C p > and 5 P > such that 

P2 (Dp(S p , s)) > C p s 
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where p,i denotes the two dimensional Lebesgue measure and 

D p (Sp, s) = {x £ D\x ■ p = h D {p) - s}. 

Then we decompose D into 

d p(&p) = {xe D\ h D (p) - 5 p <x- p< h D {p)} 

and D\D P (5 P ). The integral i T 2 e ^{x- P -h D {p)) on D\D p (5 p ) vanishes as r ^ oo. 
On D p (6p), we have 



r 2 e 2r(x- p-h D (p)) 



D P (6 P ) 



= t 2 ds e- 2rs dS 
JO JD p {5 p ,s) 

> CpT 2 [ " se- 2rs ds 
Jo 

So 1 



se 2s ds — > -C n as t — > oo. 



Proof of Lemma 4.4: The proof follows a similar scheme to [5] for Helmholtz 
equations. 
Noticing that 



f D £\F, \ 2 dx J D e|E | 2 (ix 



0(t~ 2 ) for r > 1, 



f D ^ _1 |V A Eo| 2 (ix f D p\Ho\ 2 dx 
it suffices to show 



lim 



fn\D e \ E - E o\ 2dx 



lim — = 0. 



r^oo J D p 1 \V A~Eo\ 2 dx r^oo J D p\H.Q\ 2 dx 
To estimate the numerator, we consider the boundary value problem: 

V A P = iupQ, V A Q = -iueP - in U\D, 
v A P|gn = 0, 
v A Q\dD = 0. 

or equivalently 

V A (jT x V A P) - se^V{V • eP) - cu 2 eP = eE inQ\D, 
v A P|an = 0, 
zvA(VAP)| aD = 0. 



Note that 



V • eP = 0, in Q \ D 



because V • eE = in ft \ D. 



(4.10) 



(4.11) 



ii 



Since uj is admissible, the boundary value problem (4.11) is well-posed for 
eE = e(E — Eq) which is in (Q\D). Moreover, by Proposition A. 3, one has 
P £ (H 2 (Fl\D)) 3 satisfying 

II p Wh 2 (q\d) - ^ll E IL 2 (n\o)- 
By the Sobolev embedding theorem, we have 

\P(x)-P(y)\<C\x-y\ 1 / 2 \\E\\ 



sup |P(a?)| < C||E 

Since 



l 2 (q\d)- 

xen\D 



V A (m _1 V A P) - io 2 eP = eE, 
integration by parts gives 



f _ e\B\ 2 dx = [ _ E • (V A V A P) - co 2 eP)dx 
Jn\D Jn\D 

[ _ m -1 (V A E) • (V A P) - w 2 eE ■ Pdx 
Jn\D 

[ - [ J E ■ (y A _1 V A P))dS 
J an JdD/ 



= / _V A(/i _1 VAE) -P-w 2 eE-Pdx 
Jn\D 

-([ -[ V A (^ ly ae ) - pds 

\Jan Jan J 
= - I v A (/i _1 V A E ) • PdS. 

JdD 

Expanding the RHS of the last equality at xo, one has 

f _e\B\ 2 dx = [ (P(xo)-P(x))-^A(/x" 1 VAE )fi5- f uj 2 eB ■ P(x )d 

Jn\D JdD Jd 

< cij \x-x Q \ 1 P\vA'Ho\dS + J |E |da?J||E 

< clj T\x-x Q \ l l 2 e T{x - p - hD{p)) dS + j e T( - x - p - hD ^dx\ ||E 



L 2 (n\D) 

^ l 2 (q\d) 



This yields 

/ fn\n e\E\ 2 dx\ 
[|Ef ra da, < C l MD ' ' 



< C \^U \x-x \ 1/2 e< x - p - h °Wds\ 2 + (J e^^^dxY 
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It's easy to see that the second term ( f D e T ^ x ' p ~ hD ^dx\ can be absorbed by 
the denominator J D ^ _1 |V A Eo\ 2 dx in (4.9) for large r. Therefore, it's sufficient 
to show 

lim t f \x- x \ 1/2 e r( - x - p - hD ^dS = 0. 
T ^°° JdD 

This is shown in [5], where the assumption that the Gaussian curvature of dD 
at xq is non-vanishing was used. This completes the proof of the lemma, hence 
proves the theorem. ■ 

5. Enclosing "magnetic-soft" obstacles and inclusions 

In [6], the reconstruction procedure for sound-hard obstacles also works for 
sound-soft obstacles. Inspired by this, our method also applies to enclosing a 
"magnetic-soft" obstacle. Suppose our domain O, obstacle D and all the elec- 
tromagnetic parameters in the background satisfy the same hypothesis in the 
"magnetic-hard" case, except that the fields (E,H) satisfy Maxwell's equation 
(2.1) with the boundary condition 

v A E|ao = 0. (5.1) 

Then the reconstruction scheme Theorem 4.1 applies simply by noticing the fol- 
lowing key equality. 

Lemma 5.1. Let a = 0. Assuming Maxwell's equation with the boundary condi- 
tions (5.1) and v A E|^ = v A Eo|an; we have (E, H) satisfying 

-1= /i~ 1 |VAE-VAE | 2 -cA|E-E | 2 dx+ / /i -1 ] V AE | 2 - uj 2 e\E \ 2 dx 
Jn\D Jd 

(5.2) 

where 

1 dS. 



iu (z^AE )- (u A H - v A H ) A v 
Jdn L 



Therefore, the proof essentially follows the "magnetic-hard" obstacle case, ex- 
cept that to show 

|2, 



lim / 



e|E - E \ z dx = 0, 

>Q\D 

we implement the regularity of the solution for an auxiliary boundary value prob- 
lem similar to (4.11) 

VAP = iw/iQ, V A Q = -weP - ^E in U\D, 
v A P|an = 0, 
_ uAP\ dD = 0, 

see [8] , and the Sobolev embedding theorem. 
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A. Well-posedness of a mixed boundary value problem for 

Maxwell's equations 

Let n C K 3 be a bounded domain with smooth boundary and Dcfi, Consider 
the boundary value problem of the Maxwell equation 

V A E = iujfiH V A H = -W7E in ft \ D, 
uA-E\ an = f, (A.l) 
v A (/^V A E)|aD = iujv A H| 9D = g. 

Here e and 7 are complex- valued functions in C fc (ft \ D) with positive real parts 
and uj G C. 

Theorem A.l. There is a discrete subset S o/C such that for uj not in E ; there 
exists a unique solution (E,H) G Hjj^ v (Q\D) x (Q\D) of (A.l) given any 

f G T H k j^^ 2 {d^l) and g G TH k j^^ 2 {dD). The solution satisfies 

ll E ll^(fAD) + H H H^(fAD) ^ C ^^TH^ 2 (dn) + ^\\ TH ^/2 (dD) ) (A.2) 

C > independent of f and g. 

Here we proceed to prove the theorem by modifying the variational method in 
[8], [2] and [9]. From the Maxwell's equation (A.l), the electric field E satisfies 
the second order equation 

VA(^ 1 VAE)-w 2 7E = in ft \ D, 

and 

V-7E = in ft \ D. (A.3) 

Therefore, we consider 

VA(^ 1 VAE)-s 7 - 1 V(V-7E)-o; 2 7E = in ft \ D (A.4) 

where s is a positive real number. The equation (A.3) will follow from 

sV ■ (7 _1 V(V • 7E)) + w 2 V • 7E = 

which is obtained by taking divergence of (A.4). Denoting the L 2 (ft \ D) inner 
product by (•, •) and L 2 (T) inner products by (•, -)r (where T = 9ft or dD), we 
define the bilinear form associated with the elliptic system (A.4): 

B(E, F) := V A E, V A E) + s(V • 7E, V • 7E) (A.5) 

for E, F G X where 

X = {F G (H 1 ^ \ D)f I v A F| an = 0, v ■ 7 F| aD = 0}. 

By Green's formulae, we have that B is related to the differential operator 

P = V A (^VA) - ^-^(V • 7) 

14 



by 

B(E, F) = (PE, F) - (u A (/x-V A E), F) e( „ xB) + <sV ■ 7 E, 1/ • 7F) a( ^ S) (A.6) 

for E,F G (H 1 ^ \ D)) 3 . Then for f E X' the weak formulation of the mixed 
boundary value problem 

PE = / in \ 5, 
uAE\ dn = 0, 
i/A( ) u- 1 VAE)| aD = 0. 

is: Find Eel such that 

P(E, F) = (/, F) for all F G X. 
By (A.6), this implies the natural boundary condition 

^A(/i- 1 VAE)| 9D = 0, V- 7 E| ac = 0. (A.7) 

To show the theorem, one first has for the homogeneous boundary conditions, 

Proposition A. 2. Suppose 7 and fi are complex functions in C 1 (0 \ D) with 
positive real parts, and let s be a positive real number. There is a discrete set 
S s C C such that if to is outside this set, then for any f £ X' there exists a 
unique solution E G X of 

V A (/x _1 V A E) - s 7 " 1 V(V • 7E) - u 2 7 E = / (A.8) 

satisfying 

W E Wm(n\D) ^ c \\f\\x>- 

Proof. It's sufficient to show that B is bounded and coercive on X. It's clear 
that B is bounded, 

|P(E,F)| < C'||E|| jf/1( ^ S) ||F|| H1( ^ ) . 

To show the coercivity, first we have 

|P(E,E)| > c||V AE|| L2(fA5) + s||V • E|| L2(fAS) - C , ||E|| L2( ^ S) ||E|| H1( ^ S) . 

It can be shown that there is a Poincare inequality for 1-forms in X similar to 
that in [17] 

W B \\m(n\D) - C (\\ B W 2 L2{n\D) + H V A B W%(n\D) + H V " E llL2(n\5))> EgI - 
we have 

P(E,E) > c||E||^ 1(nvS) - C\\E\\ 2 L2{n ^ y 
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To show the higher order regularity of solutions, we define for k > 2, 
X k = {F G (H k (n\D)f | vAF\ dn = V- 7 F| en = vjF\ dD = ^A(^ 1 VAF)| 9D = 0}. 
Then we have 

Proposition A. 3. Let 7 and \i be functions in C k (Q \D), k > 2, with positive 
real parts, and let s > 0. Suppose to ^ S s , then for any f G (H k ~ 2 (Q \ D)) 3 the 
equation (A. 8) has a unique solution E G X k and 

H E H// fe (fAD) - c \\f\\H k -2(n\D)- 

This can be proved by the same techniques in Section 5.9 of [17] for the Hodge 
Laplacian. 

Proof of Theorem A.l: As in [8], we take £ to be the set Si in previous 
propositions, then s can be chosen such that u ^ S s (for more details see [8]). 

To show uniqueness, suppose (E,H) G Hjy.(Q \ D) x (Q \ D) solves 

(A.l) with f = g = 0. One has 

V A (// _1 V A E) = cj 2 7E, V-7E = 0. 

It follows that Eel (by the natural boundary conditions) is a solution of (A. 8) 
with / = 0, which implies E = H = by Proposition A. 2. 

For existence, given / G TH^ /2 (dn) and g G TH^} /2 (dD) C TH^ /2 (dD), 
we can find E G iTj^. (0 \ D) with 

uAE \ dn = f, i/A(/i _1 VAE )|aD =9 
such that the extension is bounded, namely 

ll E o|lij&(n\D) - C(ll/llif=-i/2(9n) + llf llH fe - 3 / 2 (S£))) 
Suppose E G X k is a solution, given by Proposition A. 3, of (A. 8) with 

/ = -V A (yu^V A E ) + s 7 - 1 V(V • 7E0) + w 2 7 E G (H k ~ 2 (n \ D)f. 
Notice that 

ll/ll^-2(n\5) ^ C \\ E o\\ H k^\D)- 
Then E = E + E G (H k (Fl \ D)) 3 satisfies 

V A (/i _1 V A E) - S7 _1 V(V • 7E) - u 2 7 E = 0. (A.9) 

This implies V • 7E = by an earlier argument and particular choice of s. If 
we define H = A E G (H^ 1 ^ \ D)) 3 , then we have (E, H) the solution of 

Maxwell's equation. 

Applying the same argument to H, which satisfies a second order elliptic system 
by eliminating E from the original Maxwell's equation. By uniqueness, one has 
H G (H k (tt\D)) 3 . 
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The fact E <E H^. (SI \ D) is obtained by 
Div(i/AE| a(nNS) ) = -^-VAE| 9(nvS) = -^-H| a(nvS) e (H k -V 2 (d(n\D))) 3 . 
Finally, the estimate (A. 2) is derived from 

ll E ll// fc (Q\D) - H E ollH fe (Q\D) + ll E H// fc (n\D) 

< CdlEolljffc^S) + ll/lli/ fc - 2 (Q\D)) 

< C(||/|| 

H k - 1 / 2 (an) + llffllif*- 3 / 2 ^!))) (A. 10) 

< C(||/|| 
The same computation applies to H, 

H H ll#fc(fA5) ^ C '(ll/ll.ff*-3/2(afi) + |Mliffc-i/2(dD)) 

< C(ll/lliyfe-i/2(aQ) + ll5lliffc-i/2(5D))- 

Then 

||Div(f A E|a(n\S)) llHfe-V2(a(n\S)) 

< C( || Div(/) 1 1 + • H|| Hfc _i /2(aD) ) 

< C(||Div(/)|| fffc -i/ 2(an) + ||H|| ff * (nNS) ) 

< C(||Div(/) 11^-1/2(90) + ||/||jjfc-i/2( 9 n) + IMI ^-1/2(91})) 

Div Div 
The same estimate can be obtained for H. ■ 
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